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It is argued that a Hartree-Fock (HF) solution for Coulomb interacting electrons in a simple model 
of the Cu-O planes in La2Cu04, is able to predict some of the most interesting properties of this 
material, such as its insulator character and the antiferromagnetic order. Moreover, the natural 
appearance of pseudogaps in their paramagnetic and superconductor phases are also suggested by 
the discussion. These results follow after the elimination of some symmetry restrictions which are 
usually imposed on the single particle HF orbitals. One of them is the simplification of their spinor 
dependence to be of the so called a or /3 types. This constraint, seems to strongly reduce the largest 
space of orbitals corresponding to the rotational invariant HF formulation originally introduced by 
Dirac. We also remove, the demand on the HF orbitals of having a Bloch structure in the Bravais 
lattice of the crystal. This procedure allows for the consideration of HF solutions having the same 
symmetry than the known (symmetry breaking) antiferromagnetic structure of the material. It 
turns out that the most stable HF solution of the problem corresponds to an antiferromagnetic 
and insulating state (lAF), which magnetization rests on the experimentally determined direction. 
Encouragingly, the evaluated magnetic moment per cell is 0.67 /^s, a result that satisfactorily 
reproduces the experimentally measured value: 0.68 hb- Another HF solution having a slightly 
higher energy arises and corresponds to a paramagnetic state showing a pseudogap (PPG). It follows 
after only imposing the Bloch stucture on the single particle states. Finally a third paramagnetic 
but metallic solution (no gap) is also obtained by including both of the mentioned restrictions. 
The interesting result follows that this state only differs from the PPG pseudogap state in the 
form of excited empty orbitals. That is, the occupied single particle states and the HF energy at 
T = for both solutions, are identical. In general, the discussion helps to clarify the role of the 
antiferromagnetic correlations in the structure of the physics of the HTSC materials. In addition, 
these initial results, indicate a promising way for start conciliating the Mott and Slater pictures in 
the physics of the strongly correlated electron systems. 

PACS numbers: 71.10.Fd,71.15.Mb,71.27.+a,71.30.+h,74.20.-z,74.25.Ha, 74.25.Jb,74.72.-h 



I. INTRODUCTION 



The Hubbard type of models in the theory of strongly correlated electron systems are notably 
successful^'^'^'^'^'^'^'^ '^^i^-'^i-'^^'^^i-'^^'-'^^i-'^^'^^i^^'^^i^" . In particular, it is remarkable the way they reproduce the prop- 
erties of Mott insulators, such as metal-transitions oxides and copper-oxygen layered HTc compounds-'^^. However, 
the efforts for developing approaches having more basic foundations had not ceased, due to the expectation that they 
could open the way for obtaining more exact and specific results^'^'^i. In this sense, methods that are grouped into 
the so called Band Theory picture are also known as first principle calculations in the literature. They are electronic 
structure calculations that begin with the interactions among electrons or atoms in vacuum. The study of the band 
structure they predict, in principle should offer a road toward the effective and precise determination of the physical 
properties of each materiaP^'^^. Some of them are: the Configurations Interaction scheme (CI); the Local Density Ap- 
proximations method (LDA)^'^, the Local Spin Density Approximations procedure (LSDA) and Hartree-Fock method 
(HF). However, the above mentioned potentialities of those first principles approaches had been failing in describing 
many of so called strongly correlated electron systeroii. For example, the LSDA, a sophisticated generalization of 
the LDA procedure, was devised to describe local spin structures^!. However, although the method had offered sat- 
isfactory descriptions of the physical results in few materials, this success had not been universal and it also wrongly 
predicted the properties of some compounds, by example, the here considered La2Cu04. 

The motivation of the present work arose from a primary suspicion that perhaps the self-consistent Hartree Pock 
(HF) method, could had been underestimated in its possibilities for helping in the above described searches2i2iiS. In 
this sense, it can be firstly remarked that it is widespread the criteria that for obtaining behaviors such as the Mott 
insulator character, it becomes necessary the presence of short range correlations among electrons with spin quantized 
in different directions. By example, paraphrasing one type of Mott's argument for specific systems:"... two electrons 
with spin resting on contrary directions are forbidden to occupy the same Wannier orbital... ". On another hand, the 
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orthodox HF approaches does not have in consideration the correlations among electrons of different spins. Therefore 
two electrons with opposite spins do not disturb one each other and consequently both of them can occupy the same 
Wannier orbital. Clearly, the usual HF approach, seems not to be viable for investigating a system in which the Mott's 
argument is appropriate. However, the physical sources of the validity of cited Mott's statement in some systems are 
not completely clear. By example: which is the physical origin of these short range correlations assumed in it?. Even, 
the proper concept of correlations, roughly described as: "everything missing in the single-particle HF state for to be 
the real many body ground state" , makes clear how clouded their origins are. As the result of the study presented 
here, we believe that many of the so called correlation effects, can be effectively described even in the framework of the 
HF scheme, after removing certain symmetry restrictions which obstacle the finding of the best HF solutions. Such 
constraints are usually imposed on the space formed by the single particle orbitals, which are employed to construct 
the determinant like states among which the HF one shows minimal energy. By example, if after solving the HF 
problem, it occurs that the resulting self-consistent potential breaks the symmetry of the original crystalline lattice, it 
could create a gap and thus produce a Mott kind of insulating solution. This effect was originally discovered by Slater 
in Ref. 0. This symmetry breaking effect has been also more recently underlined and deepened in Ref. |23- However, 
the removal of this kind of symmetry restrictions alone had not been able to describe the insulator properties of a 
large class of materials-^ One of the central results of the present investigation, as it will be described just in what 
follows, is the identification of another important kind of symmetry restrictions that seemingly had been overlooked 
for lon g ti me. It can be cited here that a fully unrestricted formulation of the HF problem was early done by Dirac 



This work will consider the Hartree-Fock self-consistent problem as applied to a simple one band model of the 
La2Cu04SS^ but following an unusual way. In order to leave freedom to obtain paramagnetic, ferromagnetic and 
antiferromagnetic solutions in the same context, we look for single particle orbitals being non separable in their 
spacial and spinor dependence, i.e. they will have the structure 4){x,s) ^ (f>{x)'ip{s). In other words, in those states 
there is no an absolute common quantization direction for the electron spin. Thus, in each position the spin is 
quantized in a specific direction, and the equations of motion to be used will reflect this fact. Note, that to proceed in 
this way is not other thing that to apply the Dirac's unrestricted formulation of the HF procedure-^-. Wg think that the 
restriction to a and /3 types of orbitals, usually employed in HF electronic band and quantum chemistry calculations, 
prohibits from the start the prediction of possible spontaneously symmetry breaking effect o^^i^^ . Such a particular 
structure, excessively reduces the space of functions to be examined and consequently annihilates some possibilities 
to obtain exotic solutions (like that ones that are present in the strong correlation effects). We believe that in the 
context of the band theory, or more precisely, under the HF approach, it could be possible to reproduce the main 
characteristics of a wide class of Mott insulator kind of materials. The results of the simple model investigated here, 
as it will be seen, support this possibility for important compounds such as La2Cu04, which under certain doping 
levels turn out in high temperature superconductors^''. The present work heads in that direction, with eyes in also 
showing the potentialities of implementing the above considerations in self-consistent HF calculations, not only for 
describing electronic bands and the associated ground states, but also for studying atomic and molecular structures. 

The exposition proceeds as follows. In the Section |TT] we describe the details of the HF self-consistent method to 
be employed in next sections. Specifically, it will be discussed the imposition of restrictions on the space of single 
particle states in which the solutions will be searched and the possible physical consequences they could carry on. In 
Section IIIII the effective model we are going to employ will be exposed. The symmetry restrictions to be assumed 
and its corresponding tight-binding Bloch basis will be defined. Section IIVI is devoted to derive the HF equations 
associated to interacting electrons which free hamiltonian is given by the specially constructed effective tight binding 
model. The tight binding Bloch basis associated to the model is chosen in this section to have the maximal symmetry 
given by the group of translations leaving invariant the copper-oxygen plane in La2Cu04. Considering the previously 
mentioned paramagnetic solution in a generic form, the free parameters of the effective hamiltonian are adjusted to 
reproduce the form of the only half filled band in the Matheiss calculation of the band structure of La2Cu04^. After 
defined the free hamiltonian of the model, the paramagnetic and antiferromagnetic solutions are obtained. The Bloch 
basis shown in Section Hill was employed. Both solutions are also compared in this section and the properties of all the 
obtained states are analyzed in corresponding subsections. Appendix [Xl presents notations, algebraic developments, 
constants and definitions cited along the writing. In a final seccion the main conclusions of the work are reviewed and 
various further tasks for its extension are described. 



In the language of Quantum Mechanic the state of a system of N particles is described by a function depending on 
each one particle's spinor and spacial coordinates fn{xi; si, ...,xn; sn), where n represents the corresponding set of 




II. ROTATIONAL INVARIANT HARTREE-FOCK METHOD. 
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quantum numbers^^. The HF approximation consists on supposing that the above mentioned state can be expressed as 
a hnear combination of iV-products of orthonormaHzed orbitals (f>ki {xi, Si) with i = 1, N. Each one of these orbitals 
is interpreted as a single particle state, because it defines amplitude and probability distributions depending on a 
single particle coordinates. As usual, in what follows the word coordinates will mean the spacial as well as the spinor 
ones. If the particles are fermions, the previously mentioned linear combination is called Slater determinan t^^'^^i^^ . 
Let 

'H{xi,...xn) = ^T-Coixi) + ^^V{xi,Xj), (1) 

be the N-electrons system hamiltonian, including kinetic plus interaction with the environment hamiltonian TCq, 
besides Coulomb interaction among pairs of electrons V. The HF equations of motion for this system, leading the 
dynamic of the single particle states in a self-consistent way, are 

■ni s' 

-^[^ / (fx' (t)*^{x' , s')V{x, x')(j)ri{x' , s') ](t)r,^{x,s) = £^ (t)r,ix,s), 

m s' 

where ij = ki, k^^ is a label in the basis formed by the solutions. That is, each HF electron state is influenced by 
the presence of the other ones. The self-consistent hamiltonian has two components, the coulomb like type of mean 
potential which their fellows create: the direct potential and the contribution reflecting the fact that two electrons 
cannot occupy the same state: the exchange potential24. 

The HF energy of the N electrons system and the interaction energy of an electron in the t] state with the remaining 
ones, are given by 



(3) 

7] 7), 171 7), 171 

(4) 

The bracket notation definition is given in Appendix \X\ within the Subsection I A II It can be noted that the system of 
equations ([2]) is rotational invariant because it is written without imposing a spatially absolute direction for the spin 
quantization of the single electron orbitals. This rotational invariant formulation of the self-consistent HF procedure 
was firstly introduced by Dirac in Ref. [2^ . 



A. a, (3 and symmetry restrictions. 



To solve ([2]) is a complicated task because it is a system of coupled integro-differential equations. The iterative 
method is one of the most frequently employed for solving this kind of systems and it is usually complemented by 
the imposition of symmetry restrictions that simplify the space of states to be investigated. However, the use of such 
constraints could avoid the obtaining of special solutions non obeying the added symmetry conditions. Although in 
some cases they could retain the minimal energy one, the method can hide the existence of interesting excited states 
and inclusive could wrongly predict the excitation features in some cases, as it will be seen in what follows. A very 
common symmetry restriction usually employed in band theory and quantum chemistry calculations is to consider 
that single particle solutions of have spin quantized in a given direction in every point of the space^>^. That is 



(Pl.[x) ui{s) p state 



where u^i represent the Pauli spinors with spin up and down in a certain direction respectively. If the spacial functions 
and 0^ are the same, the HF calculation is called a restricted one, if they are different, the procedure is called 
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unrestricted^^. It is interesting to investigate the consequences of considering the possibly existence of non separable 
single particle states being solutions of the HF problem. A positive answer to this question can open a natural context 
for obtaining solutions exhibiting magnetic properties and to allow their comparison with paramagnetic ones. 

Another important kind of restrictions posed on the first principles band theory evaluations is the a priori impositions 
of crystal symmetries. To beforehand impose a symmetry on the assumed to be searched solution, although it could 
be well rooted in the studied system, has the risk of hiding a possible spontaneous breaking of that invariance. This 
could occurs due to the reduction of the space orbitals in which we are searching. In that case it can turns out 
that the obtained solutions will not be an absolute extremal of the energy functional, but a conditional one, due 
to the fixed symmetry constraint. For instance, let us consider the functional space formed by the allowed orbitals 
and the maximal subset of orbitals U which is invariant under a certain group of transformations T. Consider also 
the maximal subset being invariant under the group of transformations Tg which is a subgroup of T. Then, the 
set Us obtained from imposing less symmetry restrictions a priori, should contain the set U. Therefore, after to find 
the extremes of the same functional in U and Us, it could be possible to obtain different results. In this case, the 
solution in Us, in general, shall be the most stable of both. However, also it could be the case that looking for a 
solution in Us , an extremal function also pertaining to U arises as a solution. In such a situation, such a configuration 
could be found from the beginning by finding the extreme of the functional in U, that is, by imposing more symmetry 
restrictions. In terms of the HF scheme, this could mean that the states corresponding to both solutions have identical 
occupied single-particle states, but they curiously might show different sets of excited ones. Therefore, depending on 
the particular features of the material, removing a priori imposed symmetry restrictions on the set of allowed orbitals 
of the HF procedure, can predict new properties for the excited single particle states of the system. Such one could be 
for instance the gap appearance. This effect can have physical relevance after noting that at finite temperatures the 
more stable state will be preferred by the system and then the state showing a gap should be expected to be selected 
at non zero temperature. 

III. TIGHT BINDING ELECTRON MODEL: "REMOVING SYMMETRIES". 

In this section the basis of the effective band model used to describe the dynamic of the less bounded La2Cu04 
electrons will be presented. The main considerations for defining the model and the determination of its characteristic 
parameters are given. In Subsection IIII Al the simplified electronic model for the copper-oxygen planes, as well the 
main definitions in its structure are introduced. Subsections IIII Bl and IIII CI are devoted to define the symmetry 
transformation group defining the tight binding Bloch basis. 

A. Model for the Cu-O planes. 

It is known that at low temperature La2Cu04 is an antiferromagnetic-insulator— . However, in evident contradiction 
with the experiments the Linear Augmented Plane Waves (LAPW) technica^^ predicts a metal and paramagnetic 
zero temperature properties for this material. Nevertheless, such band calculation results show that the conduction 
electrons are strongly coupled to the Bravais lattice centers of the copper oxygen planes. Clearly this tight-binding 
behavior is determined by the interaction of the electrons with its surrounding effective environment. This defines 
the initial hypothesis of our model. 

The less bounded electron in the La2Cu04 molecule is the Cu^+'s not paired one. That is, at difference from O^^ 
ions, Cu ones do not have their last shell (3d) closed. Those copper 3d electrons fill the last band of La2Cu04 solid 
and in what follows they shall be referred as: the electron gas. It seems appropriate to consider those electrons as 
strongly linked to Cu02 cells and moreover, given the above mentioned arguments, with special preference for the 
Cu centers^. Thus, our Bravais lattice is going to be the squared net coincident with the array of copper sites (see 
the figure [1]) . The presence of electrons pertaining to the various fully filled bands in the material plus the nuclear 
charges, plays a double role in the model. Firstly, it will act as an effective polarizable environment, which screens the 
field created by electron charges constituting the electron gas in the half filled band. Consequently, we will introduce 
a dielectric constant e which will screen the Coulomb interaction. Secondly, as suggested by its spacial distribution 
and magnitude, the mean field created by the environment will be assumed to act as a periodic potential being 
responsible for tight-binding confining of the electron to the Cu centers. 

It is also primordial in the model to take into consideration the interaction Ff, among the electron gas and the 
"jellium" neutralizing its charges. This background will be modeled here as a gaussian distribution of positive charges 

1 y2 

/'b(y) ^ ;^exp(- — ), (6) 
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surrounding each lattice point and with characteristic radius b 
In resume, the free haniiltonian of the model takes the form 



N .2 

Ho(x) = ^-E^ + W^^(x)+F,(x), (7) 

i=l 

W.y{x) = VK^(x + R), 

where pf is the i-th electron's squared momentum operator; m is the electron mass; eq is the vacuum permitivity and 

I with and G 

moves on Bravais lattice. The versors and £^re resting on the direction defined by the lattice's nearest 
neighbours, see figure [T] a). It is known that the distance between Cu nearest neighbours is p « 3.8 A^^. We also 
consider the interaction among pairs of electrons in the form 

V{^,y) = f ^ (9) 
47reeo |x - y| 

which, as remarked above, includes a dielectric constant associated to the presence of the effective environment. 

We are seeking here for HF solutions with orbitals having a non separable spin and orbit structures. Thus, it was 
considered that the spin can show a different projection for the different Wannier wavepackets to be superposed in 
defining those orbitals. The spin for each of them will be either or a or /? type, according they are linked either to 
one or the other of the two sublattices shown in figure [T] a). Thus, the single particle eigenstates will be chosen to be 
invariant, only under the reduced group of translations which transform each of those sublattices on itself. 

It is important employ a procedure leaving independently identified the characteristics of the electron states in each 
one of the sublattices. That will allow to analyze solutions with dissimilar qualities in a same framework, that is: 
antiferromagnetic AFM, paramagnetic PM and ferromagnetic FM ones. For this purpose, let us define the points 
of the two sublattices with indices r = 1, 2, as follows 

R('^) = V2nip qi + V2n2P q2 + q^"^ (10) 
with ni and rt2 G 




where (\i and qj form the basis versors on each one of them 

a) 
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FIG. 1: The figures show: a) The point lattice associated to the Cu-O planes. For the search of the AFM properties of the 
conduction electron, and more generally for removing the symmetry restrictions, it will be helpful to separate the lattice in the 
two represented sublattices; and b) shows the corresponding base of the Cu-O planes. 
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B. Translations on the sublattices. 

The solutions we are looking for will be eigenfunctions of the operators T^cr-) belonging to the reduced discrete 
translation group which transforms a given sublattice on itself: 



= exp(i k-RM) 



(11) 



If the Bravais lattice were infinite, the Brillouin's zone (B.Z.) associated to Tj^(r) would be the shadowed one on 
figured] a). Note, that the continent square in this figure represents the B.Z. associated to the group of translations 
which leaves invariant the absolute lattice (the lattice formed by the Cu atoms in the CuO planes). However, given 
the impossibility of considering an infinite lattice for numerically solving the HF problem, then is also not allowed to 
consider its associated B.Z. as continuous. Therefore we will impose periodic boundary conditions on the 0k, z in the 
absolute lattice's boundaries xi= -L p y L p, X2= -L p y L p (see figure[2]b)). This condition determines the allowed 
set of k 



(12) 




< ni ± ^2 < 



Therefore, after recalling the discussion given in the introduction, note that we are now demanding less crystal 
symmetry on the single particle states which we are looking for, since a lower number of constraints are being 
imposed on the space of single particle states in which the solutions are searched. 






-1 -0.5 



FIG. 2: In the figures: a) The Brillouin zone associated to the infinite absolute point lattice. The vectors k label the 
eigenfunctions of the group of translations TyiM in this infinite point lattice. The grey square indicates the corresponding 
Brillouin zone (B.Z.) of the sublattices. The length of its side is ^pl-n jp. b) The net of points shows the discrete character of 
the Brillouin zone when the absolute point lattice is finite with periodic conditions fixed in its boundaries. The unit of the 
scale means -. 



C. Tight-Binding basis 

The tight binding Bloch basis which we are going to use is the following 



(13) 



^^'^^^(x,.) = J^«-^(.)^exp(zk.RW)^i,(.,(x), 



<Pr('-) (x) 



I 



, (x-RM)2 

exp( r-^ ), a<p. 



2 a2 



where N is the number of electrons in the electron gas, is the spin z projection operator, where z is the orthogonal 
direction to the copper oxygen Cu02 planes; = -1, 1, are the eigenvalue of the previously mentioned operator and 
r = 1, 2, is the label which indicates each one of the sublattices. As we are going to work on a half filling condition. 
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then N coincides with the number of cells in the crystal with fixed periodic boundary conditions Nc- Note, that due 
to the tiny overlapping among nearest neighbors approximation, the exact orthogonal character is only weakly lost 
between elements corresponding to different sublattices and having the same spin quantization. That occurs because 
nearest neighbors belong to different arrays. However the orthogonality between different elements corresponding to 
the same sublattice, as well as unity norm for every elements, is rigourously maintained. This follows because they 
are constructed as Bloch states in their corresponding sublattices. The Wannier orbitals (po{x — R*^'"') represent the 
probability amplitude of encountering one electron in the vicinity of the site R'-''-', that is on the given cell Cu02. 

Let us now describe some simplifications that will be adopted in order to solve the HF problem. The central aim 
of this work is not to make an exact study of the problem, instead we seek for approximate solutions qualitatively 
well reflecting the physical properties of La2Cu04 and other compounds. Following this principle we take for the 
Wannier orbitals the explicit form given in p3|) . Physically, this means to consider the effective potential created by 
the environment on each electron of the half filled band, is a quadratic function having a minimum on the copper 
sites and strongly confining the electrons to it. This last consideration is similar to the one made on the t-J one band 
model. 



IV. MATRIX PROBLEM AND SOLUTIONS. 



In this section the main results of this work and their discussion are presented. In Subsection IIV Al the equivalent 
matrix problem, resulting from projecting the HF system of equations ^ on the tight binding Bloch basis (fT5|) defined 
in the previous section is presented. Subsection IIV B] shows how after imposing the maximal translational symmetry 
and an a and (3 spin nature on the orbitals, the model is able to reproduce the dispersion relation of La2Cu04 
metal-paramagnetic half filled band of the precise band calculations given in Ref. 221 The solutions presented in 
subsections IIV CI and IIVD| illustrate the consequences of releasing symmetry restrictions on the space in which HF 
solutions are searched. The first of them corresponds to a Mott's insulator-antiferromagnetic ground state. That is, 
the state corresponds to an insulator even tough it has one electron per cell, precisely in the same way as La2Cu04 
behaves. The second solution obtained corresponds to a paramagnetic state showing a pseudogap and exactly the 
same HF energy and set of occupied single particle states, that the previously obtained metal-paramagnetic state^^. 
It is important to mention that all the band diagrams shown in this section are plotted on the same energy scale 
and the zero energy point coincides with the Fermi level of the antiferromagnetic ground state (lAF) presented in 
Subsection HVCl 

A. Tight-Binding representation. 

Let the searched single particle states represented in the explicitly nonseparable form 

= Y.Br:jy^^^'\^,^), (14) 

where I is the additional quantum number needed for indexing the stationary state on question, which we are going to 
define precisely further ahead. After subsituting ^ and ©, in ((21); followed by projecting the obtained result 

on the basis ip^,'"'^ and an extensive algebraic work, it is possible to arrive at the following self-consistent matrix 
problem for the coefficients appearing in the expansion (|14p : 

K + X {Gt - GZ'^ - FOl-B'^^' - e,(k) Ik-B^^', (15) 

where each of the quantities 

B"^' = Hr,L)\l (16) 

represents a vector having four components given by the four possible pairs (r, ct^). The constants 



2 2 

me a 



e^ik)^—ei{k), (17) 



8 



are dimensionless. In them, e represents the vacuum charge of the electron; h is the reduced Planck constant; a is the 
characteristic radius of the Wannier orbitals and p is the nearest neighbors separation. It is clear now that we can 
define Z = 1, 2, 3, 4, as a label indicating each of the four solutions to be obtained for every value of quasi-momentum 
k. Also, all the implicit parameters in the following 4x4 matrices are dimensionless 





II 

lPk,(t,a,) 




- 


1 1 /^dir 


;(r-,<T,) 


f^ind 
'^k — 


1 1 /^ind 
ll'-^k,(t,a^) 




Fk = 






Ik = 


P^k,(t,a^); 


ir,a-^) 1 



The set of quantities constitute the matrix representations of the periodic potential created by the mean field , 
the direct and exchange terms in the interaction potential with the neutralizing "jellium" of charges F;, defined in 
([7]) and the overlapping matrix among nearest neighbors, respectively. Each one of the four pairs, (t,az) and (r, cr^) 
defines a row and a column respectively, of the matrix in question. The explicit forms of the matrix elements are 
given in Appendix \X[ 

In this new representation the normalization condition for the HF single particle states and the HF energy of the 
system take the forms 



1 = B''''*.Ik.B'''', (19) 
<f = E ©(^^ - ^'(k))[e^(k) - |b''^''.(G*'- - Gr'').B'^-'], (20) 

where Q is the Heaviside function. 

The system (ITSt is non linear on the variables B^'J , which are the four components of each vector B''''. They can 
be interpreted as a measure of the probability amplitude of finding an electron in the state (k, 1), in the sublattice r, 
with spin z-projection az- Thinking in numerically solving the equations by the method of iterations, it is convenient 
to pre-multiply them by Ik for each k. Note that for each k four eigenvalues (/ — 1, 2, 3, 4.) will be obtained , or 
equivalently, four bands on the Z.B. (Eq. [23)1 . From Eq. ^5]) it can be observed that in the representation (fT3)) . the 
HF potentials and in general the total hamiltonian of the system, resulted as block diagonal with respect the sets of 
states indexed by the same k. This fact is a consequence of the commutation of each one of them with every element 
of the reduced group of discrete translations. That is, the group of translations which leaves invariant a sublattice. 



B. Maximally translational symmetric solutions. 



In this subsection we will search for HF solutions having their orbitals on the space of Bloch functions being 
eigenfunction of the maximal group of translations leaving invariant the absolute lattice. In other words, we demand 
the maximum possible symmetry under translations. It will follows that our model is capable of acceptably reproduce 
the profile of the conduction band dispersion calculated in Ref. [13 for the La2Cu04. We will adjust the free parameters 
of the model in order to reproduce in the best way the band dispersion reported in Ref. '221. The parameters are: 
the dielectric constant of the effective environment e; the characteristic radius of the gaussian Wannier orbitals a; the 
jumping probability between nearest sites for an electron 7 (it is fixed by the effective environment) and the radius 
in which the gaussian orbitals associated to the neutralizing "jellium" of charges decays b (see Section UlI Al) . 

In what follows the wavy hats will mean dimensionless, see Appendix [A]|A2p . Let us define the Bloch basis for the 
space of orbitals in which the solution will be searched as 

^-(x,s) = yii.-^(s)5]^?(x-R) (21) 

Fl, 

(^?(x) = exp(i Q-R) (^o(x), (22) 

where 

with n^;^, 11x2 G (23) 
and - ^ <nx^, n^^ < ^. 
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are the quasimomenta of the single particle Bloch states which are eigenfunctions of the maximal group of translations. 
Also is important to define N — L x L, and R, which are the amount of cells in the absolute lattice and the 
corresponding parametrization ([8|), respectively. The functions ipo{^) are the gaussian orbitals defined in Section 
IIII CI The searched HF orbitals will have the form 

4>Q,i{^,s) = ^sQ;'<^^'(x,s), (24) 

as expressed on the above mentioned basis. This time, the equivalent matrix problem for the "vector" i?*^"' is of 
second order for each value of (Q,l). That is, its solutions will be two component vectors. Consequently, / will take 
the values 1 or 2 now. Thus, in analogy to (|15p the new set of equations to be solved results in the form 

[E^ + X (G^^'^ - Gq'' - Fq)].b'^'' = siiQ) Iq.b'*''. (25) 

Let T and T anyone of the 2x2 matrices in ([25]) and its 4x4 equivalent on psp respectively, the relationship between 
these matrix elements is the following one 

t,r 

The direct and exchange potentials respective matrices 2x2 and 4x4 relationship becomes slightly more complicated. 
Besides also satisfying the above mentioned relation, each vector components in their "4x4" definitions must be 
removed from the sublattice label dependence and multiplied by (those new quantities are the vector components 
of the 2x2 problem). For instance, making reference to the definitions given in Appendix [XI 

-fQ,(Q,,o-,) Sa,,aAl00 + 2/oi (cOS Qip + COS Q2P) ] • 



b) epara(Q) 




FIG. 3: The figure a) shows the Brillouin zone associated to the absolute point lattice. The grey zone signals the occupied 
states in the paramagnetic metallic solution at half filling conditions. The unity of quasimomentum is ^. Figure b) shows the 
doubly degenerated bands associated to the same paramagnetic and metallic state. Note the close correspondence between 
these results and the ones obtained by Matheiss in Ref. [2^ . The zero energy level in all the band diagrams is the Fermi energy 
of the isolator and antiferromagnetic solution presented in subsection IIV CI The domain of the plot is the B.Z. of the sublattice 
shown in Fig. 2 a). 

Analogously to the ones showed in the previous section, in this representation, the normalization condition for the 
single particle states and the HF energy of the system take the forms 

1 = b'^''*.Iq.b'^'', (27) 
<f - E0(e-.-^.(Q))[^'(Q)-fB'''*-(<-<)-B^']- (28) 

The method employed for solving all the self-consistent matrix problems considered in this work was an iterative 
one which started from a particular and estimated as convenient state configuration. In the case being examined in 
this section, for beginning the iterations, we used a paramagnetic state. The figureOshows the paramagnetic, metallic 
and doubly degenerate band obtained from the iterative process on (j25p . A half filling condition has been assumed, 
that is, a state with one electron per cell is being considered. Specifically, for the case of = 20 x 20 electrons. 
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the occupied states inside the B.Z. are showed in figure [3] a) by the points inside the shadowed region. The chosen 
parameters were: e=10, which is a common value for semiconductors, a=0.25, 6=0.05 and 7=-0.03 (see Appendix [A|) . 
by following the criterium of fixing a bandwidth of 3.8 eV™. The here obtained band topologically coincides with the 
conduction band presented in Ref. j22. In both of them the Fermi level on the T-X direction is a square which vertices 
touch the middle of the B.Z. (of the CuO lattices) sides and also the maximal and minimal energies lay on coinciding 
points. Therefore, in this subsection the basis of the effective model employed here and its main set of parameters 
have been defined. 



C. Insulating and antiferromagnetic solutions. 

As we have stated before, the solution of the system of equations Eq. ([TS]) was performed by the method of 
successive iterations. The results which presented from now on were found by employing the parameter values e, a, 7 
and b, which were determined in the previous section. It is important to be noted that a, b and also that 7 must 
be of the order of the overlapping among nearest neighbors factor. It was necessary to start the iteration process from 
a particular state having an antiferromagnetic character from the beginning, in order to achieve convergence toward 
the solution presented in this subsection. In the figure [3] two sets of results following for a half filling band are shown. 
They correspond to two lattices of 20x20 and 30x30 cells. The bands are depicted on the same scale of energies. The 
difference between them is of the order of 10 dimensionless units of energy -^^= 8.3 eV. Evidently, they are bands 
corresponding to insulating states. The close similarity of both results indicates that the thermodynamical limit has 
being satisfactorily achieved for the considered sizes of the periodic system. The HE energy of this HE solution was 
the lowest among of all the ones found. In coincidence with the experimental evidence, they are states with local 
magnetic moment resting on the direction of the sublattice X12 (see figure [5]). In the next section we will show the 
difference between this HE energy of this state and the ones corresponding to the other determined HE solutions and 
shall comment on this respect. 



a)e20x2o(k) b)e30x3o(k) 




FIG. 4: Energy bands obtained for: a) A sample of 20x20 cells, Eg„p = 1.32 eV. b) A sample of 30x30 cells, Egap = 1.32 eV. 
The parameter values chosen were a = 0.25, 6 = 0.05, 7 = -0.03 and e = 10. The zero energy level is fixed in the Fermi level of 
the 20x20 system. Note that the difference in energy between the two bands is not appreciable in the employed energy scale. 
The domains of both plots is the B.Z. of the sublattice shown in Fig. 2 a) 

One important quantity which has been experimentally measured is the magnetization. It is therefore motivating 
to inspect the prediction of the obtained HE state for this magnitude. Its definition is given by the expression 

™W = E E •^kM(x, s)ct{s, s>kM(x, s'), (29) 

k',1 s,s' 

where 

cr(s, s') = (Tj;i (s, s') + a^^ (s, s') + (Tz(s, s') e^, 

and a XX = ( 1 ) ; ^^.i = ( i'i ) = ( -1 ) ^'^'^ the Pauh matrices. 

In figure [5] a) the only non vanishing component of (|29p in this solution is plotted. An interesting result is that it 
has been experimentally observed that La2Cu04 has a magnetic moment of 0.68 /is per Cu site on the CuO plane^^. 
The value obtained from evaluating the above formula for our HE state turns out to be 0.67 /is. Therefore, the 
considered here discussion satisfactorily predicts the whole antiferromagnetic structure of La2Cu04. 
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a) mi2(x) 



b) 




/ 1 
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FIG. 5: The magnetizaion vector m of the more stable HF state here determined rests in the direction 1-2. a) This figure shows 
the projection mi2 of the dimensionless magnetization, in the 1-2 direction. The magnetization unit is b) The picture 
shows a scheme of the mean magnetic moment per site in the lattice. For the shown solution its modular value is O.GT/is. 

It is also helpful to remark that the corresponding single particle states carry a more intensive antiferromagnetism 
as more closer they are from the Fermi surface. Therefore, this property offers a clear explanation of the gradual 
loss observed in the antiferromagnetic order under the doping with holes^^. In figure 6 the dependence of the angle 
^ between the magnetic moments per cell on each of the sublattice 1 and 2, shown by each one of the single particle 
Bloch states, is plotted. These components are defined as the integrals of the magnetic moment density over all 
the unit cells of the absolute lattice centered in the sublattice points. Note that the states laying just on the Fermi 
surface are perfectly antiferromagnetic ones, and that the more away from the boundaries the orbitals are, the 
less antiferromagnetic they become. Then, this HF solution indicates that after the orbitals are allowed to spatial 
dependent spin orientations, the electrons prefer to reorient their spin when traveling between contiguous lattice cells. 
This effect that can be interpreted as clean "correlation" effect when considered within a restricted HF picture.. 



FIG. 6: The single particle states exhibit a sharp antiferromagnetism in the proximities of B.Z. boundaries. In the figure the 
angle between their magnetic moment components on each of both sublattices (after divided by vr) is plotted against their 
Bloch states quasimomenta. Note that the states on the boundary have a perfect antiferromagnetism and that they become 
less antiferromagnetic as their quasimomenta move away from the boundary. The region of the plot is the B.Z. of the sublattice 
shown in Fig. 2 a) 

It also follows that the size of the zone in which the antiferromagnetism is strongest, inversely depends on the 
dielectric constant. Thus, the less is the Coulomb interaction among the electrons in the half filled electron band, 
smaller becomes the antiferromagnetic character of the single particle states and the region in which the antiferro- 
magnetism accumulate. Only the single particle states staying exactly on the Fermi surface remain having a perfect 
antiferromagnetism. In addition, the magnitude of the gap also decreases with the increasing of the dielectric constant. 



(k)/7r 



0.5 



1.0 





D. Paramagnetic solution showing a pseudogap. 



It can be observed that to remove the translational symmetry restrictions does not mean renouncing to obtain 
a paramagnetic state as solution. As a matter of fact, this solution effectively exists and its properties are quite 
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interesting. In Section |T] we had already commented about La2Cu04 and its antiferromagnetic Mott's insulator 
properties. However, at an intermediate level of doping, this material presents very special properties. After the 
breaking of the antiferromagnetic order and in a certain temperature T and doping ranges, the material transits to 
phases showing a so called pseudogap. In the La2Cu04 and other HTc materials the presence of this property has 
been observed in some regions of their paramagnetic phase PM and in the superconductor phase SC^^. 



a) 6para(k) a) PscudoGap 




FIG. 7: In the figure: a) The band structure associated to paramagnetic ground state showing a pseudogap is shown, b) A 
frontal view of the plot more clearly showing the pseudogap. The zero energy level is the Fermi energy of the lAF solution. 
The two graphics are plotted in the B. Z. of the sublattices, that is, the grey zone in figure 2 a). 



In figure [7] the band spectra corresponding to the HF paramagnetic ground state obtained from Eq. is shown. 
Note the existence of a pseudogap which reaches a maximum value of 0.1 eV w 10'^ K (equivalent to 0.012 dimensionless 
unit of energy — 8.3 eV). The parameters given in the Subsection IIV B] were employed for this evaluation. The 
pseudogap amplitude rises for smaller dielectric constants e, that is, with the reduction of the screening. It is an 
interesting result that the HF energy of this ground state is exactly coincident with the one corresponding to the 
paramagnetic and metallic state presented in Subsection IIV Bl Moreover, the lowest energy band in both solutions, 
also exactly coincide, with an upper bound error of 10"^ in dimensionless energy units (that means 10^^ in eV). Thus, 
the occupied single particle states in both solutions are identical and in consequence the momentum dependence of the 
filled energy bands also coincide. Henceforth, the difference between the two solutions refers only to the non occupied 
states. Those states can not be considered in the band model presented in Subsection IIV Bl For instance, consider 
the expressions for two states, one occupied an another empty, which are associated to the same four quasimomcntum 
value 



B 



K 






V J 



'o(x, s) = Jjj -U|(s) J2r R(-) *5Xp(i K • R^'^)) ip^(,) (x), 



X, s) = .a u^is) E..R(., (-!)'■ exp(i K • R^-^)) cp^,„ (x) 



(30) 



(31) 



where V'rC'-) represents the gaussian orbitals (po centered on R^'^-' . 

Comparing with the expression ([M)) . it may be noted that the occupied single particle state ([50]) can be expressed 
on the basis pT|) . and in fact, coincides with the state (j)J^''^ obtained from solving ([25]) . However the excited state 
(l?T|) was not allowed in the space of orbitals employed in solving (|25p . and thus its attainment is a neat consequence of 
the removing of the usually imposed symmetry under the maximal group of translations. More precisely, these excited 
states showing a pseudogap, appeared thanks to the allowed independence between the Bloch functions defined in 
both sublattices, introduced in this work. This result shows an interesting potential use that this kind of constraints 
removing procedure could introduce in the determination of new excitation properties of the HF solutions in various 
applications. 

In the following table the HF energies per particle are shown for the paramagnetic-metallic (PM), paramagnetic 
with pseudogap (PPG) and insulator-antiferromagnetic (lAF) ground states; with zero point energy assumed on the 
last one: 

It is an interesting outcome that the energy difference PM (PPG)-IAF and the Neel temperature of this kind of 
materials are both of the order of 10^ K. Thus, the results suggest the possibility of having further success in applying 
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II State 1 


IAF|| 


PMj 


PPG II 


||AE (eV)| 


0.0|| 


+0.076 1 


+0.076 II 



TABLE I: Energy differences between the various HP states. 

the approach started in this work to the description of some regions of the phase diagram of the La2 CUO4. On the 
question about the relative stability between the PM and PPG states, it can be estimated that due to the presence 
of a pseudogap, the PPG ground state should become more stable than the PM at non vanishing temperatures. This 
should be the case because for creating excitations on the PPG ground state, temperatures of hundreds of Kelvin 
degrees are needed, while excitations in the PM ground state can occur in every range of temperatures. 



Para vs. AntiFerro. 

e(k) 




FIG. 8: The figure shows in the same plot the occupied bands corresponding to the states PPG and lAF. The difference in the 
energies of the orbitals is concentrated in the boundary of the Brillouin zone. The zero energy level coincides with the Fermi 
level of the lAF solution. The domain of the plot is the B.Z. of the sublattices, given by the grey zone in figure 2 a). 

In figure [8] the PPG (PM) and lAF occupied bands are depicted in a common frame. The main difference in 
their energies corresponds to the single particle states being closer to Fermi surface. As we had noted before, the 
same behavior has the antiferromagnetic character of the single particle states of the lAF solution. Therefore, it 
can be hoped that from the solution of (j25p under doping with holes, both solutions move toward a common ground 
state without absolute magnetic order, in correspondence with the pattern shown by the phase diagrams of these 
materials^^. That is, at a critical doping level value a breaking of the antiferromagnetic order could occur and the 
system presumably transits to a phase showing a pseudogap. 



V. CONCLUSIONS 

The results of this work support the potential of the HF self-consistent method for the description of some secular 
properties exhibited by materials such as copper-oxygen compounds, specifically La2Cu04. Those are properties 
usually associated to strong correlation effects, that could be explained by first principle calculations, after modified 
to incorporate procedures to search for spontaneous symmetry breaking solutions and general ways of exploring the 
spin structure of the HF obitals. Note that single particle orbitals are not only employed in the HF method. They 
are also essential ingredients of more general schemes as the various types of density functional methods. Thus, the 
kind of reasoning employed in this work seems to be easily easily implemented in such discussions. In the present 
work we explore those paths in a simplified manner. In order to avoid the intrinsic complexities of La2Cu04 material, 
it was necessary to utilize a simple model which was sufficiently flexible for reproducing the dispersion profile of the 
only half filled band of La2Cu04 reported in Ref. [2^. After imposing the maximal symmetry under translations a 
paramagnetic and metallic ground state (PM) was obtained which dispersion properties topologically also coincide 
with the results given in Ref. Then, the free parameters of effective model were fixed from the requirement 

of reproducing the bandwidth of the half filled band obtained in Ref. [221. Employing those parameters and from 
removing some symmetry restrictions, solutions were obtained which show interesting properties. In agreement with 
the experiment, the insulating-antiferromagnetic (lAF) solution turns to be the most stable among all HF states 
found. Some of its properties are enumerated below: 
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1. The gap magnitude diminish with the increasing of the screening constant e. 

2. The antiferromagnetic structure the HF orbitals increases when the states approach the Fermi level. In addition, 
the size of the outlying region in which antiferromagnetism persists depends on the screening created by the 
effective environment e. That is, by increasing screening, the size of the antiferromagnetic region reduces. Thus, 
the idea arises that after doping with holes (that is, solving for HF solutions not at half filling condition as it is 
done here) the antiferromagnetic zone, which is precisely concentrated near the Fermi level could be annihilated, 
producing in this way a phase transition to a non magnetically ordered ground state. This possibility indicates 
a way to describe the normal state properties of the HTc superconductors through a simple HF study. 

3. The magnetic moments per cell which are evaluated rest on the right direction respect to the lattice and its 
modular value of 0.67 /is satisfactorily coincides with the experimentally measured result 0.68 hb- This outcome 
support the adjust made of the parameters defining the effective model employed. 

The other HF solution presented in this work corresponds to a paramagnetic state showing a pseudogap (PPG). 
Few properties of this state are described in what follows: 

1. The pseudogap magnitude diminishes with the increasing of the screening constant e. 

2. Just as we hoped from the beginning of the work, at zero temperature the states PM and PPG resulted identical. 
The difference between them is only given by the excitations of the system. Thus, the removal of some symmetry 
restrictions, defines in this case, new properties for unoccupied single particle states. It seems feasible that in 
other materials, it could be possible to obtain a gap instead of a pseudogap even in the absence of magnetic 
order. Such an outcome could show the ability of a properly formulated HF description in describing all the 
kinds of Mott insulators, being or not magnetically ordered. 

3. The amplitude of the obtained pseudogap reaches a value of the order of lO^K. Thus, it is possible that the 
system at finite temperature should prefer the PPG state over the PM one. That is an interesting signal, 
knowing that La2-a;M2,Cu04 presents pseudogap in the PM phase. 

4. In the same way as it happens for the antiferromagnetic character of the lAF ground state, the difference between 
the one particle energies of lAF and the PPG (PM) states, is greater for states being closer to the Fermi surface. 
That is, it happens for the electrons with more energy and consequently the first ones in disappear under doping 
with holes. Therefore, as it was mentioned before, these results suggest the possibility to describe a transition 
lAF^PPG under doping with holes, within the here studied effective model. 

One methodological conclusion of this work seems to be worth to be mention. It corresponds to the fact that the 
results fully clarifies that the solutions of a general HF problem non necessarily should turn to be a set of single 
particle orbitals, all having an a or /3 spin structure at all points of the space. This is directly shown by the particular 
example of the lAF solution in which the HF single particle states have a neatly non separable character in their spin 
and orbital dependence. This possibility, although being simple, and perhaps even expected to be correct after to 
direct the attention to it, looks to be relevant for the description of magnetic and strongly correlated systems. 

The above enumerated conclusions motivated new objectives to be considered in extending the work. Our general 
interests are twofold: in one side to achieve a more complete and precise investigation about the here detected 
unexplored potentialities of the first principle calculations. In another direction, to look for the possibility of describing 
the HTc superconductivity in the framework of the here studied simple model. Some of the more specific issues of 
future searches in connection those general objectives are the following ones: 

• To generalize the discussion done in this work to introduce the doping with holes as a new parameter. Then, 
it will be possible to investigate the effects of the doping on the here determined ground states. Of particular 
relevance in this sense appears the lAF state. 

• To compute the zero temperature electron Green function of the system, and use it to evaluate the effective 
polarization of the La2Cu04 in the obtained states. 

• With the polarization results at hand, attempt to solve the Bethe-Salpeter equation for two holes in the lAF 
ground state (and also in the PPG one), to find whether it is possible or not to define the existence of p reformed 
Cooper Pairs in the model. The possibility for its existence was suggested by the results of Ref. [30, in which 
it was argued that a strong 2D-screening of the Coulomb interaction is created by a half filled band of tight 
binding electron. 

• To continue the application of the ideas advanced here in help clarifying the widely known debate between the 
Mott and Slater pictures about the electronic structure of solids. 
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• To apply a rotational invariant HF calculation for obtaining states of molecular or atomic systems having 
incomplete shells. We guess, that this application could help to reduce the value of the correlation energies 
in those systems. This idea is suggested by the results of this work which show that the HF energy can be 
optimized in many of the currently considered physical systems, due to the difficulties in determining the best 
among all possible self-consistent solutions. 
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APPENDIX A: MATRIX ELEMENTS. 
1. Brackets notation 

The bracket terms in ([3]) represent the following integrals: 

(Al) 

(miTiob) = ^ d(2;,s) Ti-aix) (/>p(x,s), 

s 

{m,n\V\o,p) EE ^ (fx(fx'(f)'^{x,s)(f)'^{x',s')V{x,x')(f>o{x',s')(f)p{x,s), 



where m, n, o and p, denote anyone of the possible quantum number arrays. 

2. Dimensionless definitions. 

Before writing the expressions for the operators matrix elements, we begin by showing some other useful definitions, 
as for instance the employed dimensionless ones: 

Unit of distance, (A2) 
Dimensionless characteristic length, (A3) 

Dimensionless lattice points position, (A4) 

V Dimensionless Coulomb potential. (A5) 

From those definitions, the used reference conventions can be inferred. That is, the distances are expressed in units 
of the Cu nearest neighbors separation p = 3.6 A; and the energies and potential interactions, in units of the quantity 
that for the employed parameters is equivalent to 8.3 eV. 

3. Other definitions and properties. 

We had already defined along the paper the Wannier orbitals ^nir), as normalized gaussian functions, centered in 
R^*"^ and with characteristic parameter a. In order to simplify notation when the jellium in ([7]) is going to be used. 



p = 


1 




a 


a = 






P 




R 


R = 






P 






V = 


ft2 
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is useful to define: 

^S(x) = ^exp(-i^). (A6) 
Let us use (|A6P for defining the jellium's potential 

f,(x) ^J2[d'y ^R(y)v^(x-y)^^(y). (A7) 

R 

In obtaining the matrix elements of the jellium potential (|A7[) . the following notation is useful: 

(rW,rV|R',R^*^) ^ I d'xd'y ^*^,^,{^)^Uy)V{^-y)v'k{y)^Ri^^{^), (A8) 

where labels r and t move on the sublattices independently. In the context of the infinite lattice problem is easy to 
see the following property 

(R(''^R''lF|R^R^*^> = (r'"^ -R,o''|F|o^R^*^ -R). 

(A9) 

which can be proved by making a pair of changes of variables on the second order integral in the right hand side of 
(|A8|) . For finite systems the prove is a little more complicated. Firstly, the integrals extend over the region occupied 
by the lattice, and for this region the integration does not remain invariant under the translations that must be done. 
It becomes necessary to periodically extend the Coulomb interaction beyond the boundaries. That is, modified it to 
the form 

y„(x,y) = y -4 r-^> (AlO) 

where ni and n2 ~ -1, 0, 1, and L is the length of the side of the squared region which the system occupies. 
However when the system is sufficiently large, the error done by no extending it periodically is not important; and it 
vanishes in thermodynamic limit. 

Equally useful in defining the direct and exchange potential matrix elements, is the following notation 



(RW,R(*')|F|R(*"),rW) EE y'd2^d2y(^*^,.,(x)^^,„,(y)y(x-y)^j,„»,(y)¥'RU)(x), (All) 

where t' and t" also move on both sublattices independently. In the same manner they fulfill the following property 
(rM,R(*')|F|R(*"\r(*)) = (R('') -R(*"),R(*'' -R(*")|F|0,R(*) -R(*"'), (A12) 

if we define 

p^(t',t") ^ R(t') _R(t"). (A13) 

Then it follows 

^ - \r(2), ift'^t". ^^^^> 

Let us consider the definitions of R^^^ and R'-^-' given in pU|) . and recall that the only Wannier orbitals which have 
non vanishing overlapping are those centered on the same site or the ones centered in nearest neighbors (the closest 
neighbors belong to different sublattice). Then for fixed R'^'"-' and R*^*-*, the only no vanishing among all the quantities 
in the left hand side of (|A12[) are 

(A15) 

(R(^) -R(*"\o|1/|0,R(*) -R(*")> for t'=t". 
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and 

(A16) 

(R('^) -R(*"\pjy|0,RW -R(*")) for t' T^t", 
in which The four quantities Pi are defined as 



P e^i, 


if 


*=1, 




if 




P ^X2 1 


if 






if 





That is, they move over the neighbors which are closest to the site on the origin. 

The employed procedure reduces the number of integrals appearing in (|A15|) , (|A16p , inclusive those corresponding 
to the right hand side in (jA9P and anyone of those appearing when the matrix elements of the periodic potential , 
or the projection of the tight binding Bloch basis between any two elements, are searched. By example 

(R^'-^IM^^IR^*)) = j (fx <p*^,^,W^ ip^w, (A18) 
(rM|rW) ^ j fx^*^,^, (A19) 

which respectively fulfill the following properties 

(R(")|VF^|R(*^) = (R("'*'|VF-,|0), (A20) 
(R^'-^IR^*)) = (R^^'^^IO), (A21) 

given the periodicity of W-y in the absolute sublattice. 

In the following subsection it is frequently used the symbol 5r^t+i in which t + 1 is not the usual sum of 1, but he 
transformation of a given sublattice in another 

[2 ift^l, 

t + l = I ' A22 

I 1 if t=2. ^ ' 



4. Matrix Elements. 

Making use of the previously definitions given in this Appendix and after performing an extensive algebraic work, 
the desired matrix elements are computed. Below we start presenting them: 

^k,{t,a,),(rM,) = ^a,,a, Wqo St,r + '2rf{cOS kip + COS k2p) 5t,r+l , (A23) 
4,(t.a,),(r-,<T,) = da,,aAl00 St.,r + '2Ioi{cOskip + COSk2p) St,r+l], (A24) 
-Pk,(t,a,),(r,<T,) = Sa,, cr A^OO St,r + '2Foi{cOS kip + COS k2p) St.r+l], (A25) 
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where Wqo = 0; represents a change in the zero point energy; and 7 is a free parameter describing our lack of knowledge 
about the periodic potential. The other appearing parameters are defined as 

100 - (0|0) (A26) 
= 1, 

101 = (0|Pi) (A27) 

= e~^, 

Foo = ^X^(R,0''|y|0^R), (A28) 

R. 

= ^E(^ + Pi 'O'll^lO^R), (A29) 

R, 

where N is the number of electrons in the electron gas and, as we had already defined, the symbol ~ means dimen- 
sionless. 

The matrix elements of the direct potential are 

(A30) 

k',1 

^0*'*"^ - I E (r''*"\o|v^|o,r'''"^, (A31) 



where 



Z^^''"^ ^ ^Y. H cos(k.p,) (p,; + R^*''"\o|y|0,R^*'*"\ (A32) 



^ 1^ ^ cos(k.p, + k'.p,.)(P. + R^*''"\p,l^|0,R^*''"V (A33) 



For simplicity, in expression (|A30P the Einstein summation convention for the indices {t' , t" , a'^ , a'^) is employed. 
Further, the matrix elements of the exchange potential are 



(A34) 



^ind _ r nkM* ^i^M' ,t,r) „kM , nk',!* <-,(k.k',t,r+l) „kM 

k',1 



where 



pk'.l* r.(k.k',t+l.r) pk'.l ok',r r.(k,k'.t,r) pk'.l 



^(k,k',',") = 1 ^ cos[(k-k')- R(*'*")]x(R^*'''\o|y|0,R^*''"\ (A35) 



A^ , 



^(k,k',',") = 1^ ^ cos[k.p, + (k-kr R(*'*")]x(p, + R^*'*"\o|y|0,R^'''"\ (A36) 
St''^''"''"^ = ^Y. E cos[k.(p,+p^O + (k-kO- R(*'*'')]x (p, + R^*''''\p^.|F|0,R^*''"^ (A37) 
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5. Reducing the order of some integrals. 

Any one of the fourth fold integrals presented in previous sections, can be partially integrated in quadratures, in 
such a way they final calculation is reduced to numerically evaluate first order integrals. Thus 

exp[-2S^] ^ r [(^x, + |)sin</.-(i?,, -i)cos<^]2. 



{n + p^,p^\V\0,R) = ^^=^x d<t> eM- ' ^' \~2 —} (A38) 

^ r r [(^..2 - i)sin0- (^^, + i)cOS0]2 

X Erfcj 1, 

where Erfc is the complement error function. Similarly 

,~ ~ h,~, h ~, expf-;^] /■^'^ ^ [(^a,, + i) sine/)- ^a,, cos 0]^^ 

{^ + V.,0^\V\0\B.) = X I del. exp{- ^^ } (A39) 



^ — -"^ c 



^ 17 f r [^^2 sin ^ - (i?^, + i) cos 

X Ertci (i+c2)a2 }, 



where ( = 



— a 
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